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Abstract 

We investigate QCD effects in the production of a single slepton at hadron colliders in the 
Minimal Supersymmetric Standard Model without R-parity. We calculate the total cross sections 
and the transverse momentum distributions at next-to-leading order in QCD. The NLO corrections 
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CN , enhance the total cross sections and decrease the dependence of the total cross sections on the 



' factorization and renormalization scales. For the differential cross sections, we resum all order 

> 

soft gluon effects to give reliable predictions for the transverse momentum distributions. We also 
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ly-^ \ slightly different at the Tevatron and are in good agreement at the LHC. Our results can be useful 



compare two approaches to the non-perturbative parametrization and found that the results are 



to the simulation of the events and to the future collider experiments. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) is one of the most popular new 
physics models beyond the Standard Model (SM). In the MSSM, R-parity conservation is 
imposed in order to keep the proton stable. R-parity is a discrete symmetry, which is defined 
to be Rp = (_i)3{-B-i)+2S^ where B, L and S are baryon number, lepton number, and spin, 
respectively. All SM particles have -Rp = 1, and all superpartners have Rp = —1. In 
consequence, the superpartners can only be produced in pair, and the lightest superpartner 
(LSP) is stable. However, proton decay may be avoided by assuming B ot L conservation, 
but not both. In these cases, R-parity can be violated. For a recent review of the R-parity 
violating (fip) extension of the MSSM, see Ref. 

The most general form for the fip part of the superpotential is 

Wj^^ = ^.UH, + ]^\,,L,L,El + KjkUQM + (1) 

where Li{Qi) and Ei{Ui,Di) are, respectively, the left-handed lepton (quark) SU{2) doublet 
and right-handed lepton (quark) SU{2) singlet chiral superfields, and Hi 2 are the Higgs 
chiral superfields. The indices i,j, k denote generations and the superscript c denotes charge 
conjugation. As stated before, we require that the lepton number violating part and the 
baryon number violating part should not exist simultaneously. For the purpose of this paper, 
we will assume that only the lepton number violating couplings X[j^. are non-zero. A recent 
summary of the bounds on the couplings in Eq. (Q) can be found in Ref. Q. 

The /2p-MSSM has the remarkable feature that sparticles need not be produced in pair. 
Therefore, a definitive signal of R-parity violation will be the observation of single sparticle 
production. The terms in Eq. (0) will result in various resonant sparticle production pro- 
cesses. The term involving Xijk leads to resonant sneutrino production at e^e~ colliders 
while the X'-ji^ term leads to resonant squark production in hadron-hadron collisions jj, 5|- 
The term with X'--, can induce both resonant squark production at ep colliders & and reso- 

nnn 

nant slepton production at hadron-hadron colhders 12|, |8|| . We will consider the last case 
in the following. 

The resonant production of a single slepton can lead to interesting phenomenology at 
hadron colliders. A charged slepton can decay into a neutralino and a charged lepton, and 
the neutrahno can subsequently decay into a charged lepton and two jets via A' couphngs. 
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Due to the Majorana nature of the neutrahno, the two leptons can have either opposite or 
same charges. The case of two leptons of the same charges is more interesting due to the 
absence of large SM background. For a detailed analysis of the signal and the background, 
see Ref. For a sneutrino, the decay products can be a chargino and a charged lepton. 
With the subsequent decay of the chargino, the final state could contain three charged 
leptons which is also a clean signature. 

The analysis in the above is all based on the tree-level cross sections. It is well known 
that in order to reduce the scale dependence of the cross section and get more precision 
predictions, one must include the QCD radiative corrections. The next-to-leading order 
QCD corrections to the total cross sections have been calculated in Ref. P]. However, they 
used rather old parton distribution functions in the numerical calculations, and their results 
have large dependence on the PDF sets used. Moreover, they did not consider the kinematic 
distribution of the events, which is very important in designing the strategy of discovery. In 
this paper, we will calculate the transverse momentum distribution of the slepton at next- 
to- leading order in QCD, and also calculate all order soft gluon resummation effects to give 
reasonable predictions. In our numerical evaluations, we will use the new sets of PDFs, and 
compare our predictions for the total cross sections with those of Ref. j^. 

The paper is organized as follows. Section |Tll describes the production processes at 
hadron/parton level and shows the leading order (LO) cross sections. Section UTTl presents 
the analytical expression of the next-to-leading order (NLO) corrections to the cross sec- 
tions. Section HVl gives the analytical and numerical results for the total cross sections and 
Section El gives the transverse momentum distribution. And finally is the summary. 

II. THE PROCESSES AND THE CROSS SECTIONS 

We consider the process A{pi) + B{p2) L{q) + X, where A and B are the incoming 
hadrons (proton-antiproton for Tevatron or proton-proton for LHC) with momenta pi and 
P2-, L stands for charged slepton / (or sneutrino u) with momentum q. At parton level, 
four kinds of subprocesses can be induced by the X[ji^ couplings at tree level: dk + Uj — >• Ij, 
dj + dk ^ Vi, and their charge conjugated processes. 

In the QCD improved parton model, the hadronic differential cross section can be factor- 
ized into the convolution of the partonic differential cross sections with appropriate parton 
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distribution functions (PDFs): 

where ct"^ is the cross section for the partonic subprocess a{pi) + P{p2) — > L{q) + X. pi 
and p2 are the momenta of the incoming partons a, (5. The momenta fractions Xi and X2 
are defined by Xi = pi/pi, {i = 1,2). fp/H{x,fif) is the parton distribution function which 
describes the probabihty of finding a parton p with momentum fraction x inside the hadron 
H at factorization scale /if. The sum is over all possible initial partons which contribute. 
For any momentum q, it can be decomposed as 

Q = {q^,q~,QT), qT = QT, (3) 

where q^ = q^ ± q^, and the rapidity y is defined as 

1 0+ 

y = n^^-- (4) 

z q 

We work in the center-of-mass frame of the colliding hadrons, in which 

pi = ( v^, 0, Ot) , P2 = (0, v^, Ot) , (5) 

where s — (pi +^2)^ is the center-of-mass energy squared. 

The partonic cross section a'^^ can be calculated order by order in perturbation theory. 
The leading order contribution is the Born level process a{pi) + P{p2) — > L{q). Here a and 
/? represent only quarks and anti-quarks. The matrix element is simply 

Mb = X'ijMP2)PLu(pi) (6) 

where — (1 — 7^)/2, and we have ignored the mixing between the left-handed and 
right-handed sleptons since they are almost degenerate. The couphng X'-j^ depends on the 
generations of a, j3 and L. For simplicity, we will make the subscript implicit below and 
simply write A'. The one particle phase space can be written as 

= {2'K)dqldy5{ql)5{pt - q^)S{p^ - q'). (7) 
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After convoluted with the PDFs, the differential cross section with respect to the transverse 
momentum and rapidity of the slepton is 

dcJR TT 



d^ = —^\^BfS{q^)U/A{x^^,fif)fp/B{xifif) 

= KQl)fa/A{xl l^f)U/B{xl fif), (8) 

where x^ = y/rey and = y^e^^ with r = m^/s, and m is the mass of the slepton. We 
have made the summation over the initial states a, (3 implicit. 



III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO QCD corrections consist of the following contributions: the exchange of virtual 
gluons and the corresponding renormalization counterterms, the real gluon emission sub- 
processes, the gluon initiated subprocesses, and the contributions of AltareUi-Parisi (A-P) 
splitting functions. In the following, we will calculate these contributions seperately. We 
use dimensional regularization (DREG) in c? = 4 — 2e dimensions to regulate all divergences 
and adopt MS renomalization and factorization scheme to remove the ultraviolet (UV) and 
infrared (IR) (including soft and coUinear) divergences. 



A. Virtual gluon exchange 



Evaluating the relevant Feynman diagrams of virtual corrections, we obtain the amplitude 



Mv = MB7z^CF{^7ryr{l + e) 

ZTT 



1 I, l^r 1 

— ^--In— -- 



2 271^ ■ 
+ ^ - 1 



(9) 



where Cp — A/3 and Hr is the renormalization scale. In order to remove the UV divergences 
in Aiy, a renormalization procedure must be carried out. We define the renormalization 
constants as 



A'o = ZxX'i^l = (1 + SZx)X'i^% 

1 + ^SZ^l) tpL+ M + ^SZ^R ] ^R, 



(10) 



(11) 
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where Aq and ipo are the bare couphng and the bare quark wave function, respectively. In 
the MS scheme, these renormahzation constants are fixed to be 



5ZA = -^C^(47r)T(l + e)-. 

Thus the running of the couphng constant is governed by 

A'(/xo) 



(12) 
(13) 



1 . 2"^^ 1 ^ 
1 + — Cplog — 

47r /x^ 



2 ■ 



(14) 



After adding the counter term, the UV divergences in My are cancelled, but the IR 
divergent terms still persist. The corresponding differential cross section is 



da 



V 



das 



dqj.dy dqj.dy % 



Ci.(47r)T(l + e) 



1 



1/3 , /i^ 
+ ln-^ 



e V2 



2 m2 3 



(15) 



B. Real gluon emission 



We now consider the contribution from the emission of one gluon Q;(pi) + /9(p2) — L{q) + 
g{k). We define the Mandelstam variables as 

s = (pi +P2)^ 

t = ipi — qf = — pfq~ = — ^ s{q^ + m^)e~^, 
u — {p2 — (if — — P2 — — ^ s{q^ + m^)e^, 

S = (pi +^2)^ = X1X2S, 

t — {pi — qf — Xi{t — m^) + m^, 

u ^ {P2 — of ^ X2{u — w?') ^- rn^ . (16) 

The squared matrix element (with spin and color summed and averaged) can be expressed 
as 
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y 



ill 



-TTasy'^nfxiX2S^ {(I - e) 
9 q^ 



1 + 



r 



X1X2 



+ 2e- 



X1X2 



(17) 
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The two particle phase space in d dimensions is 



87rr(l -e) V?T 
Thus the differential cross section can be written as 



--l^(r~^qTdy 5{s + t + u- m^) 

dgrdySis + i+u-m^). (18) 



2(27r) 

1 1 /in 



dq^dy 



J dxida;2F(xi, X2)S{xiX2S + Xi{t — m^) + X2{u — m^) + m^) 
^ :F{x^,x*), (19) 



where 



+ + 

Xi — T ^ ^ X1X2 — T 



and 



xt^ey^T + q^/s, xt = e-y^T + qys, (20) 



X — <! (l-e) 



1 + 



-V 



XIX2 J 



T 



+ 2e }. (21) 

X1X2 



C. Gluon splitting subprocesses 

In addition to the real gluon emission subprocess, there are also contributions from the 
gluon initiated processes a{pi) + g{p2) L{q) + /3{k) and gipi) + Pip2) L{q) + a{k). 
Defining the Mandelstam variables as before, for the first subprocesses we can write the the 
squared matrix elements as 



I . .<Z9|2 1 ^ x/2..26 (l-^)(^^ + ^^)+2e^ffl^ 

l-^G I = o^"^-^ 7: TTJ 

—[1 — e)st 

1 2 2e 1 [ {X2xt - T)[{x2xt - rf + T^] 2e (3:20;+ -r)V 
^-naA A^..i.2.^| 13^ ix.X2r 



(22) 



7 



while the squared matrix elements for the second subprocess can be obtained from the above 
one by the substitution t il, xi X2, xf X2 , namely 



1 ,2 26(1 +^ ) +2etm" 



3 * -il-e)su 



3 \ {x\x-:^'^ 1 - e 



Thus the differential cross section is 
ddG 



where 

F (xi, X2) = — V — 



X 



falA{Xx,lif)fglB{x2,lif)\M^if + fglA{Xx,llf)f(ilB{x2,\^i)\M^^f . (25) 



(23) 



J dxidx2F' {xi, X2)S{xiX2S + Xi{t — m^) + X2{u — m^) + m^) 
^ :i^'(xi,x;), (24) 



D. Splitting function contributions 

The contributions from splitting functions are essential to cancel the coUinear divergences 
in the cross sections. They arise from the renormalization of the PDFs at NLO. In the MS 
scheme, the renormalized PDFs can be expressed as 

/„/^(x,/.,) = /./^(x) + ^(-i)^r(l + 6)f^^ f^'^P^p{z)U/H{x/zl (26) 



where Pap{z) are the A-P splitting functions, which are given by 



1 



2^^^-^^ + (r^, 



The resulting contributions to the differential cross section are 



(27) 
(28) 



^ = ^A-,(,2)l^r(l + e)^ 
aq^ay 12s e 27r \ 



X [(Po/),/^(a;?,/x/)/^/B(4,/x/) +/c,/a(x?,/x/)(Po/)^/b(4,/x/)] , (29) 
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where 



(30) 



IV. TOTAL CROSS SECTIONS 



The total cross sections can be obtained by integrating out qt and y from the differential 
cross sections given above. The leading order result is 

"'^ dx 



7T f dx 

(^LO = Y^^' y —fa/Aix,fif)fp/Bir/x,fif). 

The next-to-leading order result can be written as 



(31) 



ctnlo = / dxidx2 



+ a''^fa/A{xi,fXf)fg/B{x2, /i/) + a^^^ fg/A{xi, fXf) ff3/B{x2, ^^f) 



(32) 



where 



-X 



9s 

+ 2(1 + z- 



In 

2^ I ln(l - z 



1 + z 



1 - Z 



[l + Z^ 



\nz 

1 - z 



a^9 = 



48 s' 



■A 



12 



In !!L + In ) (^2 ^ _ + z){lz - 3) 



(33) 

(34) 



Tfi T 

Here z = —r- = , and the function with a subscript is a distribution, which is 

S X1X2 

defined as 

>1 rl 



(izg{z)f{z) 



dz{g{z)-g{l))f{z) 



(35) 



for an arbitrary function g{z). Our expressions for the total cross sections are the same as 
those obtained in Ref. 

In the numerical evaluation of the total cross sections, we use the updated version of the 
CTEQ and MRST PDF, namely CTEQ6.1 and MRST2004 (there is no update for GRV98 
PDF). The renormalization and factorization scale are taken to the mass of the final state 
slepton, yUj. = /ij = m, unless otherwise specified. We will set the initial quarks to only the 
light flavours u, d and s. Due to the severe constraints on the products of two A's, we will 
consider the contributions from one single A^j^i. at a time. The most stringent constraints 



on the couplings relevant here come from the rare semileptonic decay K ttuu, which 
gives X[ji^ < 0.012mj^/(100GeV) for k = 1,2. We can assume that the first two generations 
of sfermions have nearly degenerate masses, so we can take the down squark mass as the 
slepton mass approximately. In the following, unless otherwise specified, we will take the 
tree level coupling X'-jj^ = 0.01, which can satisfy the constraints for all the slepton masses 
under consideration. The cross sections for other values of A' can be obtained easily by a 
multiplicative factor. 

For convenience, we define three cross sections as following: 

(Tloi • LO partonic cross section convoluted with LO PDFs; 

(Tlo2 '■ LO partonic cross section convoluted with NLO (MS) PDFs; 

o"nlo • NLO partonic cross section convoluted with NLO (MS) PDFs, 

and correspondingly two K factors: 

ctnlo cnlo 
Ki = , K2 = . 

C"L01 <^L02 

As the above definitions, K2 measures only the size of the NLO QCD corrections to the 
cross sections, while Ki accounts for the effects of changing parton distribution functions 
additionally. 

A. Sneutrino production 

We will consider sneutrino production first. The possible partonic initial states are dd, 
ds, sd and ss. 

Fig. n shows the cross sections for sneutrino production at the Tevatron as functions 
of the sneutrino mass. The left and right graphs correspond to CTEQ6.1 and MRST2004 
PDFs, respectively. The solid, dashed and dotted curves represents to ctnlo; ctloi and crL02, 
respectively. For pp collision, since we take the A's to be equal, the cross sections from 
sd initial states are obviously equal to those from ds initial states. We find that the NLO 
QCD corrections generally enhance the total cross sections. The cross sections decrease 
monotonically with the increasing of the sneutrino mass. From Fig. one can see that, for 
a sneutrino of mass m = 200 GeV, if we take A' = 0.02, the leading order (LOl) and the 
next-leading-order cross sections from dd channel are 274 fb and 400 fb, respectively. The 
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cross sections from ds channel are 58 fb and 83 fb, while the ones from ss channel are 11 fb 
and 16 fb, respectively. The cross sections for the LHC case are shown in Fig. |21 For pp 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

mj>(TeV) mp(TeV) 



FIG. 1: The cross sections for sneutrino production from various initial states at the Tevatron. 
The cross sections from sd initial states is obviously equal to those from ds initial states for pp 
collision if the corresponding A's are equal. 

collider, the cross sections from ds and sd initial states are no longer equal even with equal 
couplings. For a 200 GeV sneutrino and A' = 0.02, the cross sections from dd channel can 
reach 3.17 pb (LOl) and 4.2 pb (NLO), respectively. The cross sections from other channels 
are relatively smaller, but still remarkable as the cases at the Tevatron. Our results are 
consistent with those from Ref. 

To quantify the enhancement of the total cross sections by the NLO QCD corrections, 
we show the K factors in the following. First, we plot K2 as functions of the sneutrino mass 
in Fig. El The curves are quite similar in spite of different PDF sets and different initial 
states, which is just a reflection of the flavour blindness of QCD. The solid and dashed curves 
correspond to CTEQ and MRST PDFs, respectively. For each PDF set, the curves from 
top to bottom represents dd, ds, sd and ss initiated processes, respectively, while for the 
Tevatron case sd is omitted. The ordering of the magnitudes of the K factors of different 
initial states is expected, which is due to the different weights of the contributions from 
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FIG. 2: The cross sections for sneutrino production at the LHC. For pp colhsion ds initiated process 
has larger cross section than sd initiated one if the couphngs are the same. 

gluon initiated subprocesses. The corrections increase monotonically with the increasing of 
sneutrino mass, and can reach 70 percent at Tevatron and 30 percent at LHC for a 1 TeV 
sneutrino, respectively. 

The case of Ki is a bit more comphcated due to the different order of PDFs involved. In 
Fig. 0] we plot Ki for dd initial states as functions of the sneutrino mass. We find that there 
is still certain discrepancy between the results of CTEQ and MRST parton distributions, 
especially at large iriy region. This is mainly due to the difference of the leading order 
results at region of large momentum fraction, while the next-to-leading order results of the 
two PDF sets are quite close to each other, as shown in Fig. ^ This confirm the necessity of 
calculating the NLO corrections. For the su 
Ki for ss initial states. From Fig. 5 in Ref. 



jprocesses involving strange quark, we only plot 
Q|] , one can see that there was large discrepancy 



between the results of CTEQ5 and MRST98 PDF sets. With the updated PDFs, we find 
that the discrepancy has been significantly reduced. As shown in Fig. the two curves 
agree each other quite well, and are more close to the CTEQ5 results given in Ref. j|9|. 

We now investigate the dependence of the cross sections on the renormalization and 
factorization scales. We will choose the renormalization /z^ and the factorization scale /x/ to 
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FIG. 3: K2 for sneutrino production at the Tevatron (left graph) and at the LHC (right graph). 
The sohd and dashed curves correspond to CTEQ and MRST PDFs, respectively. For each PDF 
set, the curves from top to bottom represents dd, ds, sd and ss initiated processes, respectively, 
while for the Tevatron case sd is omitted. 

be equal, = fJ'f = /^yUo, where /^o = fnc, is the mass of the sneutrino. In Fig. IHl we plot 
the ratios of the cross sections R = cr(/t)/cr(l) as functions of k for m,> = 200 GeV. We only 
use CTEQ PDFs, since the parametrization of the parton distribution is irrelevant for the 
main conclusion here. We show here the results of dd initial states, and the results of other 
initial states are similar. From the figure one can see that the dependence on the scales is 
significantly reduced from LO to NLO, as we expected. 



B. Charged slepton production 

We now turn to the production of charged sleptons. The possible initial states now 
become ud and us for sleptons with positive charge, du and su for ones with negative 
charge, respectively. For simplicity, we will only consider contributions from ud and du initial 
states, and only use the CTEQ6.1 PDF. The situation is much similar to the production of 
sneutrinos, so we do not discuss any details here and just show the cross sections in Fig. [7| 
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FIG. 4: Ki for sneutrino production from dd initial states at the Tevatron (left graph) and at the 
LHC (right graph). 
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FIG. 5: Ki for sneutrino production from ss initial states at the Tevatron (left graph) and at the 
LHC (right graph). 
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FIG. 6: The renormalization and factorization scale dependence of the leading order and next-to- 
leading order cross sections at Tevatron (left graph) and at LHC (right graph). 

and the K factors in Fig. |H1 respectively. 



V. TRANSVERSE MOMENTUM DISTRIBUTION 

In this section we investigate the transverse momentum distribution of the slepton. At 
Born level the slepton is kept at zero qx due to momentum conservation and the distribution 
is proportional to 5'^{qT)- Thus the leading order distribution at non-zero qt belongs to 
(9(as), where momentum conservation is retained by the additional parton emitted. The 
distribution can be obtained by integrate out y from the differential cross section given in 
Section imi However, this fixed order result is only valid when qt is not too small compared 
with the mass of the slepton m. If qx -C m, the corresponding parton emitted would 
be either soft or coUinear to one of the initial partons. In consequence, large logarithms 
like ln(m^/g|,) will appear and will dominate the cross section for sufficiently small qx- In 
general, there should be double logarithms for each gluon attached to the initial quarks due 
to the overlap of soft region and collinear region. As a result, the perturbative expansion 
would be controlled by ln^(m^/g|.) rather than a^. The convergence of the perturbation 
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FIG. 7: The cross sections for charged slepton production at the Tevatron (left) and the LHC 
(right) . 
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FIG. 8: The K factor {Ki) for charged slepton production at the Tevatron (left) and the LHC 
(right). 
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series will be spoiled if ln^(m^/g|,) approaches unity. 

The problem arising at small qt can also be seen from the fact that the fixed order 
cross section is singular as qx — > 0. For later use, we define the asymptotic part of the 
differential cross section at small qt to be the terms which are at least as singular as l/q^ 
when qT 0. At the asymptotic expression can be isolated from the expressions 

obtained in section ITTTl 

d(T 



(asym) = \(yo — \\fa/A{x\,iif)fp/B{xl,iif) (^2(7^ In ^ - 3C^^ 
z, Tx qrp \^ \ qT J 



(36) 



where ctq = 77;- A'^. 

12s 



A. Resummation of large logarithms 

In order to make use of the perturbation theory with the existence of large logarithms 
at each order, one must reorganize the perturbative expansion to resum the large terms. It 
was first shown by Dokshitzer, Diakonov and Troyan (DDT) 10] that in the double leading 
logarithm approximation (DLLA), terms like a"ln^"'~^(m^/g|.)/g|i can be resummed into a 
Sudakov form factor. However, this result relies on the assumption that the emitted gluons 
are both soft and coUinear and their transverse momenta are strongly ordered. With the 
over-constrained phase space, the resulting distribution is over-suppressed at small g-r- This 
implies that the subleading logarithms are also important and need to be resummed too. 
The subleading logarithms correspond to the phase space configurations in which some of 
the emitted gluons are not so soft or coUinear. In this case, the transverse momentum con- 
servation must be iraposed, which is implemented in the 6-space formalism introduced by 
Parisi and Petronzio Collins and Soper jl^ improved the 6-space formalism to resum all 
the terms like a^U (r = 0, ...,2n — 1) based on the renormalization group equation (RGE) 
method, where L represents the large logarithms. In this framework, Collins, Soper and 
Sterman (CSS) "l^ derived a resummation formula for the transverse momentum distribu- 
tions of the vector bosons produced in Drell-Yan processes. This formalism, often refered 
to as the CSS formalism, has also been applied to many other processes. 
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In the CSS formalism, the differential cross section we are considering can be written as 

do" da 



resum) + Y{qT, m, Xi,X2), (37) 



dq^dy dq^dy 

where the resummed part can be expressed as an inverse Fourier transformation 
^ (resum) = X] 2^°27r / ^^^^^P ' ^ap{h, m, x?, 4) 

= V -^0 / bdbMbqT)W^f,{b, m, X?, x°), (38) 

where b is the impact parameter conjugating to qt, Jq is zero order Bessel function of the 
first kind, and 

„o ^0 



W^pib, m, x^, x^) = U/Aix^i, C:,/b)fp/B{xl Cs/b) 



X exp 



dfi' 



In ^^A(a,(/i)) + B{as{ji)) 



(39) 



Here Ci{i = 1,2, 3) are constants of order 1 which are by convention [13] chosen to be 

Ci = C3 = 2e-''^ = bo, C2 = 1, (40) 

and / is the convolution of the PDFs and the coefficient functions C 

/•I 

fa/h{x,^) = ^ —Ca-y{z,a,{fi))f^/h{x,fi), (41) 

and the coefficients A, B and C can be expanded to series in 

00 

A(a.) = (42) 



n=l 

00 



5(a.) = 5^i?W(^)", (43) 

n=l 

00 

C^./3(^,a.) = 5^Ci;)(.)(^)", (44) 

n=0 

and they can be calculated order by order in perturbative theory. The lowest order coeffi- 
cients can be extracted from the asymptotic expression above 

A^"^ = Cf = \, B^'^ = -\C, = -2, (45) 
C;^i{z)=6^,6{l-z). (46) 
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With these coefficients, we can actually sum up all terms like a^L^""^ and OgL^""^. 

The another term in Eq. (jHTf) . the Y term, is the remaining contributions which are not 
resummed. Since it contains no large logarithms, it can be reliably calculated in perturbation 
theory. In fact, according to the definition of the asymptotic part mentioned above, one can 
immediately see that 



Y = J 2^ (pert) — ^ (asym). (47) 



do" da 
dq^dy dq^dy ' 

However, the resummed part is still not able to be calculated perturbatively. The reason 
is that in Eq. (jHHj) . the integral over the impact parameter b extends to infinity, while the 
integrand involves the strong coupling constant and the PDFs at scale bo/b, where they 
are not well defined if b is large enough so that bo/b enters non-perturbative region. 

Many prescriptions have been developed regarding this problem, and we will show two 
approaches here. Collins, Soper and Sterman, in their original paper [l^, suggested that 
one can use a cut-off 6max and regard the effects from b > femax as non-perturbative input. 
Practically, they replacing W{b) in Eq. (jHHj) by 

W{b) = iy(6.)FNp(6), (48) 

where 



K = , (49) 
and -Fnp(6) parameterizes the non-perturbative effects. Since 6* never exceeds Vax, W^(fe*) 



can be calculated perturbatively, and the theoretical uncertainty main 
tion Fnp. Recently, Landry, Brock, Nadolsky and Yuan (BLNY) 



y relies on the func- 
proposed the form 



NP 



exp I -b^ + In ^ + ^71^73 ln(100a;?x°) j | . (50) 

They take 6max = 0.5GeV~^, Qo = l-6GeV and the parameters giii = 1,2,3) are fitted to 
the available Drell-Yan data, which are given by 

(71 = 0.21, ^2 = 0.68, g, = -0.60. (51) 

One of the largest disadvantages of the 6=,, prescription is that it alters the W function 
in the perturbative region b < &max- Regarding this, Qiu and Zhang (QZ) [l^ proposed 
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another prescription so that W{b) is kept unchanged for b < b^g^^- Namely, they chose 



, , , W{b) b < 6^ax, , , 

W{b) = { (52) 
W{b^,,)F^p{b) b>b. 



where Fnp takes the form 



FMb) = exp |-ln— 1^ [g, {{b'r - (Cx)") + 92{b' - Cx)] - Mb' - &Lx)| • (53) 

The parameters gi and a are fixed by the continuity of the first and second derivative of W 
at 6 = 6max, while g2 and g2 are determined from experiments. In the numerical evaluation, 
we will use both above methods and compare their results. 



B. Numerical results 

In the numerical evaluation of the transverse momentum distribution, we take the renoma- 
lization and factorization scale for the fixed order expressions to be = yU/ = a/w^ + 
For simplicity, we only show the results for dd initial states as an example. The results for 
other initial states are similar. We use the CTEQ PDFs throughout this section. 

According to the previous analysis, the differential cross section with respect to and 
y can be formally written as 

da da . . da , , da 

(total) = — (resum) + —5-— (pert) — (asym). (54) 
dq^dy^ ' dq\dy^ ' d(^^dy^ ' dq\dy^ ' ^ ' 

At small qTi the asymptotic part and the perturbative part cancel each other, and the 
resummed part dominates the distribution. At large gy, the difference between the resummed 
part and the asymptotic part belongs to higher orders in a^, so the perturbative predictions 
are recovered. Since the perturbative part and the asymptotic part can be reliably calculated 
within perturbation theory, the only ambiguity comes from the resummed part due to the 
non-perturbative issue. So we shall compare the results of the different approaches to the 
non-perturbative parametrization for the reliable predictions. 

With the two 6-space prescriptions at hand (BLNY and QZ), we first compare the inte- 
grand of the Bessel transformation bW{b\ We plot this function for a sneutrino with mass 
m = 200 GeV and rapidity ?/ = in Fig. El The solid curves represent the BLNY results, 
while the dashed and the dotted curves correspond to the the QZ prescription with two 
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choices of g2- 92 = (no power corrections) and g2log{m^b'^g^^/bl) = O.SGeV^ (fitted power 
corrections), respectively. The parameter ^2 is always taken to be zero. From Fig. El we 
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FIG. 9: Integrand of the 6-integral (without the Bessel function) as a function of b for m = 200 GeV 
and y = at the Tevatron (left) and the LHC (right). 

can see that at the Tevatron, the BLNY parametrization notably changes the shape of the 
function for b < 6max, and shows a sharper decrease than the one of QZ. As a result, the 
resummed differential cross sections for the two approaches will differ at small qt (at larger 
qt the impact of large b region gets smaller due to the Bessel function). On the other hand, 
the QZ parametrization itself is insensitive to the value of g2. At the LHC, as shown in 
the right graph of Fig. IHl the difference is quite small, and the resummed cross sections are 
expected to agree to each other. Furthermore, the shape of the function is "narrower" at 
the LHC than the one at the Tevatron, i.e., the contributions from large b region to the 
integral at the LHC are smaller than the ones at the Tevatron. So we conclude that the 
predictions at the LHC is more reliable than the ones at the Tevatron. In Fig. EH we plot 
the corresponding differential cross sections according to the two approaches. As discussed 
above, the results are slightly different at Tevatron and agree quite well at LHC. In the 
following calculations, we will adopt the QZ prescription and take g2 = 0. 

In Fig. ^2 we plot the various parts of the differential cross section in Eq. for rapidity 
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gr(GeV) 5T(GeV) 

FIG. 10: Comparison of the resummed differential cross sections for m = 200 GeV and y = at 

the Tevatron (left) and the LHC (right). 

y = at the Tevatron and the LHC. The mass of the sneutrino is taken to be 200 GeV. The 
soUd, dashed, dotted and dot-dashed curves correspond to the total result, the resummed 
part, the perturbative part and the asymptotic part, respectively. Quantitatively, the per- 
turbative and the asymptotic cross sections agree very well at small transverse momentum. 
On the other hand, the resummed and the asymptotic part are not cancelled completely at 
high Qt due to the higher order effects included in the resummed one, so that the total one 
and the perturbative one will differ at large qt- This can be considered as the theoretical 
uncertainties. In principle, one can return to the perturbative result for qt > q'™*, where 
is arbitrarily chosen in the intermediate qt region. However, in order to make the transition 
smooth, one must introduce some kinds of matching procedure which could also lead to 
uncertainties. In our work, we will use the resummed cross section from small qt to large 

Having settled down the technical issues, we now turn to show the predictions for the 
transverse momentum distributions of the slepton. For the y integral, the kinematic range 
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FIG. 11: The distribution with respect to the transverse momentum and the rapidity of the 
sneutrino from dd channel for central rapidity y = at the Tevatron (left graph) and the LHC 
(right graph). The mass of the sneutrino is taken to be 200 GeV. 



should be 



\y\ < arccosh 



1 + r 



(55) 



In Fig. ^1 we plot the transverse momentum distribution of the sneutrino from the dd 
channel at the Tevatron and the LHC separately. The mass of the sneutrino is taken to be 
200 GeV, 400 GeV, 600 GeV, respectively. The peaks of the distribution appear at about 
3 GeV at the Tevatron and about 5 GeV at the LHC. The differential cross sections decrease 
sharply with the increase of qt, which indicates that most events will happen in the relatively 
low Qt region, where the resummation effects are essential. For example, about 80 percent 
of events lie in the qt < 20 GeV region at the Tevatron, while at the LHC the range is 
qt < 40 GeV. Therefore, the precise theoretical predictions for the distribution in the small 
qt region are very important, which can help the selection of the the experimental cuts for 
better background rejection. 
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FIG. 12: The transverse momentum distribution of the sneutrino from dd channel at the Tevatron 
(left graph) and the LHC (right graph). The mass of the sneutrino is taken to be 200 GeV, 
400 GeV, 600 GeV, as labeled in the figure. 

VI. SUMMARY 



In summary, we have calculated the next-to-leading order total cross section and trans- 
verse momentum distribution of a single slepton in the MSSM without R-parity, including 
all-order soft gluon resummation effects. For the total cross section, our results are consis- 
te„t With *e .suits of Ref. g We update the. ,.,„e*a. .esu.ts w.th the updated version 
of the parton distribution functions and find that the discrepancy between the different sets 
of PDFs has been decreased compared with the previous version. Also, we find that the 
factorization and renormalization scale dependece in next-to-leading order is much smaller 
than the leading-order results. For the transverse momentum distribution, we resummed 
the large logarithms at small transverse momentum. Combined with the fixed order calcula- 
tions, we give consistent predictions for both small qt and large qx- We also compared two 
approaches to the non-perturbative parametrization and found that the results are slightly 
different at the Tevatron and are in good agreement at the LHC. Our result can be useful 
to the simulation of the events and to the future collider experiments. 
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